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INTRODUCTION 

The  principal  advantages  of  laminar  over  turbulent  boundary  layer  flow  for 
operational  vehicles  are  reduced  frictional  drag  and  diminished  self  noise. 

Such  prospects  of  increased  performance  have  tempted  designers  and  engineers 
over  the  years.  However,  due  to  the  uncertainty  of  transition  and  the 
•vucial  differences  in  designs  appropriate  to  either  laminar  or  turbulent 
boundary  layer  flows,  laminar  flow  and  transition  have,  for  the  most  part,  been 
designed  around  rather  than  exploited.  The  strongest  deterrents  to  practical 
application  of  laminar  flow  have  been:  (1)  The  inability  to  achieve  transition 

length  Reynolds  numbers  of  practical  interest,  and  (2)  a profound  lack  of 
understanding  of  the  transition  process  itself,  including  prediction  of  its 
onset. 

Important  break-throughs  in  recent  years  have,  however,  attracted  more 

attention  to  and  rekindled  interest  in  these  problems.  Specifically,  considerable 

progress  has  been  achieved  in  extending  transition  to  higher  Reynolds  numbers 

utilizing  the  stabilizing  mechanisms  of  pressure  gradient^  (vehicle  shape), 

2 3 4 

surface  heating  * and  suction  . New,  elaborate  computer  codes  which  model 

5 

the  present  transition  prediction  theories  are  now  operational  and  available 
for  general  use.  These  are  being  corroborated  as  well  as  improved  by  an 
increasing  experimental  data  base^’^.  In  short,  recent  advances  in  understanding 
of  the  transition  process  (though  still  incomplete)  and  of  techniques  for  the 
maintenance  of  laminar  flow  to  significantly  higher  Reynolds  numbers  have 
brought  about  a resurgence  of  hope  for  the  practical  and  reliable  design  of 
laminar  flow  vehicles. 

Ultimately,  however,  before  a laminar  flow  vehicle  can  be  deemed  suitable 
for  application,  a reckoning  of  the  desirable  features  of  laminar  flow  against 
reliability  and  effective  costs  will  have  to  be  made.  The  risk  and  consequences 
of  early  transition  must  be  carefully  assessed.  Effective  costs,  on  the  other 
hand,  would  take  the  form  of  system  complexity,  of  limits  imposted  on  the  range 
of  operating  conditions,  and,  of  course,  of  dollars  and  cents.  The  subject 
of  the  present  report,  while  fundamental  and  technical  in  nature,  provides 
valuable  information  to  eventually  enable  this  necessary  reckoning  or  trade-off. 
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As  previously  mentioned,  three  parameters  have  been  shown  to  favorably 
affect  the  stability  and  delay  transition  of  a laminar  boundary  layer;  viz, 
favorable  pressure  gradient,  surface  temperature  (over  ambient)  and  surface 
suction  velocity.^  Increasing  any  of  these  parameters  individually  will  aid, 
to  a limit,  in  stabilizing  a laminar  boundary  layer.  When  used  in  combination, 
the  effects  undoubtedly  become  interrelated  and  more  complex.  Still, 
theoretical  studies  have  led  to  the  prediction  of  large  extensions  of 
transition  length  Reynolds  numbers  when  heat  is  added  to  the  surface  of  a 
vehicle  forebody,  where  pressure  is  decreasing  in  the  streamwise  direction, 
that  is,  in  a region  of  favorable  pressure  gradient.  These  results  have  been 
confirmed  (at  the  least)  qualitatively  in  recent  laboratory  experiments  and 

O 

in  the  Autonetics  flow  tube. 

The  focus  of  the  present  study  is  primarily  on  the  afterbody,  rather  than 
the  forebody  region.  Here  the  pressure  gradients  are  strongly  adverse  and  the 
effects  of  suction  must  be  called  into  play  in  order  to  stabilize  a laminar 
boundary  layer.  Before  describing  the  intent  and  scope  of  the  present  study, 
it  would  be  well  to  first  complete  the  description  of  its  motivation  by  means 
of  the  following  sample  calculations:  The  vehicle  configuration  shown  in 

Figure  1 (designated  as  AR-3,  one  of  the  two  utilized  as  a test  bed  in  the 
present  study)  has  a length  of  55  feet,  a maximum  diameter  of  11  feet  and  a 

■3 

displacement  volume  of  2446  ft  . Given  a power  plant  with  2000  SHP,  drag 
estimates  for  a fully  turbulent,  unheated  hull  boundary  layer  indicates  a 
vehicle  speed  of  35  knots.  Now  if  the  forebody  is  heated  with  a portion  of 
the  waste  heat  from  the  propulsion  system  such  that  laminar  flow  is  maintained 
to  the  streamwise  location  of  minimum  pressure  (approximately  72  percent  of 
length),  drag  is  reduced  and  vehicle  speed  is  increased  to  55  knots.  At 
this  condition,  the  last  28  percent  of  the  vehicle,  the  afterbody  region, 
which  is  supporting  a turbulent  boundary  layer,  contributes  about  90%  of  the 
remaining  drag. 
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If  the  afterbody  were  laminarized,  say,  by  the  addition  of  the  remaining 
waste  heat,  of  suction,  or  of  some  combination  of  the  two,  to  extend  the 
laminar  boundary  layer  to  90  percent  of  length,  the  drag  falls  to  12  percent 
of  that  of  the  fully  turbulent  vehicle  and  the  vehicle  speed  increases  to 
nearly  75  knots.  From  another  viewpoint,  the  55-knot  speed  for  the  heated 
laminar  vehicle  may  be  obtained  with  a 33  percent  reduction  in  power  for  the 
case  of  a laminarized  afterbody.  This  is  indeed  an  appealing  and  - within  the 
frame  work  of  the  available  stability  theory  - a completely  feasible  prospect. 


The  stabilizing  effects  of  suction  are  considerably  stronger  than  those 
g 

of  surface  heating  , particularly  in  a region  of  adverse  pressure  gradient. 
However,  by  clever  and  efficient  utilization  of  propulsion  system  waste  heat, 
surface  heating  appears  to  be  a less  complicated  and  more  economical  technique 
as  compared  to  suction  for  practical  applications.  With  a view  towards  the 
aforementioned  necessary  reconciliation  between  performance  gains  due  to  laminar 
flow  and  corresponding  costs,  it  is  of  interest  to  attempt  to  maximize  the 
utilization  of  heat  and  thereby  reduce  the  amount  of  suction  necessary  to 
laminarize  afterbody  boundary  layers. 


The  problem  addressed  in  the  present  study  is  that  of  the  effectiveness 
of  surface  heating,  of  suction  and  of  combinations  of  the  two  in  the  afterbody 
or  adverse  pressure  region  of  laminar  vehicles.  The  study  was  purely  analytical, 
although  support  and  guidance  was  taken  from  relevant  experimental  studies 
whenever  possible.  In  this  regard,  the  Autonetics  Flow  Tube  at  Colorado  State 
University  has  proven  to  be  a unique,  flexible,  relatively  inexpensive  facility 
for  high  Reynolds  number  laboratory  experiments.  The  present  problem  has  been 
idealized  in  the  sense  that  no  notice  was  taken  of  the  probable  engineering 
difficulties  associated  with  neatly  injecting  heat  into  the  vehicle  surface  and 
with  operating  micro-pore  suction  surfaces  on  the  presence  of  naturally  occurring 
(clogging)  contaminants.  These  further  aspects  are,  in  fact,  presently, 
simultaneously  under  study  at  Autonetics. 
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Two  hull  shapes,  both  of  fineness  ratio  5 to  1 , were  selected  on  the  basis 
of  their  widely  differing  pressure  distributions  to  serve  as  test  beds  for  the 
present  analytical  study.  These  are  shown  in  Figure  2.  The  external  or  potential 
velocity  distributions  corresponding  to  a vehicle  speed  that  yields  a length 
Reynolds  number  at  minimum  pressure  location  of  at  least  250  x 10^  (in  water  at 
40°F)  are  shown  in  Figure  3.  Calculations  were  performed  utilizing 
axisymmetric,  finite-difference  boundary  layer  computer  codes  to  determine  such 
quantities  as: 

1.  The  forebody  temperature  distributions  and  heating  requirements  for 
laminar  flow  stabilization  to  the  point  of  minimum  pressure. 

2.  The  limits  of  extension  of  laminar  flow  into  the  afterbody's 
adverse  pressure  gradient  with  surface  heating  alone. 

3.  The  continuous  suction  distributions  required  to  maintain  laminar 
flow  over  the  entire  afterbody. 

4.  The  intermittent  heat  plus  suction  distributions  which  provide 
similar  stabilization. 

5.  Stability  characteristics  for  the  various  cases  tested,  including 
critical  disturbance  frequencies  and  corresponding  amplitude  growth 
histories. 

fl 

Moreover,  total  drag  was  estimated  for  the  two  hull  shapes  under  the  various 
conditions  studied.  Wake  drag,  drag  equivalent  of  suction  and  drag  equivalent 
of  heating  have  been  accounted  for.  A total  drag  summary  serves  as  a basis  for 
comparison  of  effectiveness  of  the  different  laminar  flow  stabilization  schemes,  not 
only  among  themselves  but  also  between  the  two  hull  shapes. 

Some  discussion  of  the  theoretical  details  which  describe  the  effects  of 
pressure  gradient,  surface  heating  and  suction  on  boundary  layer  flows  is 
presented  in  Section  I.  The  step  by  step  computational  procedures  as  well  as 
an  overview  of  the  range  of  cases  investigated  is  also  given.  The  calculated 
results  are  presented  and  described  in  Section  II.  Section  III  is  dedicated 
to  the  hull  drag  analysis  while  discussions  and  conclusions  are  contained  in 
Section  IV. 
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I.  ANALYTICAL  DETAILS  AND  PROCEDURES 
A.  REVIEW  OF  ANALYTICAL  DETAILS 

The  effects  of  pressure  gradient  and  wall  suction  on  both  laminar  separation 
and  transition  have  been  under  study  for  many  years  now.  In  fact,  Prandtl 
demonstrated  the  delay  of  separation  by  suction  in  experiments  carried  out  over 
70  years  ago4.  In  the  case  of  suction,  delayed  separation  most  likely  corresponds 
to  delayed  transition  because  the  effects  of  suction  on  the  boundary  layer  are 
powerful  and  are  felt  across  the  entire  layer.  Within  the  last  decade  or  so,  the 
additional  technique  of  surface  heating  for  a water  boundary  layer  has  been  added 
to  the  list  of  boundary  layer  control  mechanisms.  The  details  of  the  analyses  and 
\ experiments  are  of  little  Interest  here.  However,  a brief  listing  of  the  major 

findings  along  with  a few  graphs  which  summarize  the  relevant  parametric 
relationships  for  stabilization  of  laminar  flow  would  perhaps  be  an  appropriate 
way  to  set  the  st^ge  for  the  present  study. 

ff 

1 . Pressure  Gradient 

It  has  been  well  established  experimentally  as  well  as  analytically  that 
, a favorable  pressure  gradient  stabilizes  and  an  adverse  gradient  destabilizes  a 

laminar  flow.**  Mechanistically,  a favorable  gradient  of  pressure  tends  to  thin 
the  boundary  layer  and  thereby  lower  its  thickness  Reynolds  number,  while  an 
adverse  gradient  destabilizes  the  flow  through  the  formation  of  an  inflection 
point  within  the  velocity  profile.  Recent  theoretical  (approximate)  results  by 

• Q 

King  suggest  that  increased  stability  produced  by  "reasonable"  changes  in  a 
favorable  pressure  gradient  are  smaller  than  those  produced  by  corresponding 
• changes  in  either  surface  heating  or  suction.  The  problem  is  extremely  complex, 

however,  in  the  sense  that  the  relative  strengths  of  the  mechanisms  depend  on 
the  actual  level  of  each.  In  other  words,  the  magnitude  of  effects  of  heating, 


for  example,  depend  on  the  magnitude  of  the  pressure  gradient  and  the  level  of 
suction.  Still,  the  above  statement  is  believed  to  be  true  over  the  parametric 
ranges  of  interest.  The  effects  of  pressure  gradient  on  boundary  layer  stability 
compete  with  the  other  two  in  a way  that  is  little  understood  - particularly 
for  the  case  of  an  adverse  gradient,  as  identified  and  addressed  in  the  present 
study.  The  effects  of  pressure  gradient,  suction  and  surface  heating  on 
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separation  of  laminar  boundary  layers,  however,  have  been  previously 
investigated.  The  results  for  separation  are  both  interesting  and 
illuminating  for  the  phenomenon  of  stability/transition  and  will  be  discussed 
in  later  sections  of  this  report. 

2.  Suction 

The  stabilizing  influence  of  suction  occurs  through  two  separate  but 
related  effects:  Thinning  of  the  boundary  layer  and  alterations  of  the  profile. 

Both  are  of  a globa1  nature  in  that  they  are  felt  throughout  the  boundary  layer. 
Suction,  therefore,  has  a powerful  effect  on  stability;  stability  characteristics 
are  intimately  dependent  upon  both  the  boundary  layer  velocity  and  curvature 
profiles.  It  is  apparently  possible,  at  least  intuitively,  to  stabilize  a 
laminar  boundary  layer  against  any  adverse  pressure  gradient  with  sufficient 
suction. 

3.  Surface  Heating 

The  effects  of  surface  heating  on  laminar  stabilization  are  less 
well-known.  The  linear  theory  stability  calculations  of  Wazzan  and  Smith  which 
incorporate  surface  heating  have  become  quite  popular  recently  and  their 
agreement  with  experiments  using  heated  models  is  indeed  impressive.  The  effects  o 
surface  heating  for  water  (Prandtl  number  of  8)  have  been  shown  to  manifest 
primarily  through  viscosity  variations  and  to  be  confined  to  a thin  layer  near 
the  wall.  Slight  alterations  of  the  velocity  profile  near  the  wall  are  shown 
by  the  calculations  of  Wazzan  and  Smith  to  have  significant  effects  on 
transition  in  both  favorable  and  adverse  pressure  gradients.  The  localized 
influence  of  heating  as  compared  to  the  global  influence  of  suction  on 
transition  suggest  that  heating  would  be  strongly  over  shadowed  by  suction. 

This  issue  is,  of  course,  central  to  the  present  study.  Few  known  calculations 
and  very  little  experimental  data  are  available  for  the  case  of  an  adverse 
pressure  gradient,  excluding  the  output  of  the  present  study. 
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4.  Boundary  Layer  Shape  Factor 

Experience  has  demonstrated  that  the  effects  of  pressure  gradient, 
surface  heating  and  suction  can  be  correlated  to  a critical  Reynolds  number 
based  on  displacement  thickness  as  a function  of  gross  boundary  layer 
parameters.  The  shape  factor,  H,  defined  as  the  ratio  of  the  displacement 
to  the  momentum  thicknesses,  has  been  used  quite  successfully  in  this 
capacity.  Figure  4 from  Reference  6 shows  the  relationship  between  critical 
Reynolds  number  based  on  displacement  thickness  and  shape  factor,  H,  for 
flows  involving  all  three  mechanisms.  To  a first  approximation,  the  critical 
Reynolds  number  is  a function  of  H alone.  This  is  not,  in  general,  strictly 
true.  More  accurate,  detailed  calculations  for  Falkner-Skan  wedge  flow 
profiles  carried  out  recently  by  Gazley  and  Wazzan^  indicate  that  the 
correlation  requires  two  parameters,  one  of  which  may  be  H.  Furthermore, 
the  history  or  evaluation  of  the  boundary  layer  - as  evinced  by  the 
streamwise  distribution  of  H - may  also  be  of  importance.  In  the  same  study, 
they  also  constructed  a first  order  correlation  of  a transition  Reynolds 
number  (also  based  on  displacement  thickness)  to  H for  the  same  wedge  flows 
utilizing  the  "e^"  criterion  to  indicate  transition  (see  Figure  5).  Their 
results  are  in  good  agreement  with  calculations  performed  by  Wazzan  for 
heated,  laminar  flow  vehicles  and  with  preliminary  experimental  results 
obtained  during  tests  with  models  of  the  same  heated  vehicles.  The  studies 
by  Gazley  and  Wazzan  are  continuing  and  quite  promising,  in  the  presence  of 
a growing  experimental  data  base  for  comparisons  to  theoretical  predictions. 

The  immediate  utility  of  such  correlations  lies  in  their  use  for  vehicle 
1 2 

design  studies  and  for  indicating  the  probability  of  boundary  layer  transition 
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without  the  necessity  of  costly,  detailed  stability  calculations.  It  was 
in  the  latter  role  that  good  use  was  made  of  the  established  correlation  between 
critical  Reynolds  number  and  shape  factor  (Figure  4)  in  the  present  study.  The 
shape  factor  is  relatively  easy  and  inexpensive  to  monitor  for  various  vehicle 
flow  conditions,  say,  by  means  of  finite  difference  boundary  layer  calculations. 
Only  those  flow  conf igurations  that  survive  scrutiny  of  their  H distribution 
need  be  investigated  for  boundary  layer  stability  by  means  of  the  relatively 
elaborate,  linear  stability  analysis  computer  code  (the  TAPS  code). 
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The  question  of  whether  a given  laminar  boundary  layer  is  stable  can  not 
be  answered  by  studying  the  shape  factor  alone,  as  presently  understood. 

A universal  transition  Reynolds  number  correlation  has  yet  to  be  established. 
Exceeding  the  critical  Reynolds  number  for  a given  H is  a necessary,  but  not 
sufficient  condition  for  transition  - hence,  the  need  ultimately  for  detailed 
stability  calculations.  More  will  be  said  about  these  factors  in  subsequent 
sections. 

B.  PROCEDURES  FOLLOWED  DURING  STUDY 

Two  hull  shapes  were  selected  for  the  present  study  the  first,  the  so-called 
AR-3,  generated  by  the  Parsons-Goodson  routine^  and  the  second,  a Reichardt  body, 
as  illustrated  in  Figure  2.  Both  have  a fineness  ratio  of  5 to  1 and  an 
assumed  length  of  55  feet.  Their  corresponding  potential  velocity  distributions 
are  shown  in  Figure  3.  Using  the  relationship 

Cp  = 1 - (U/Uco)2 

The  curves  in  Figure  3 may  be  rescaled  to  potential  pressure  coefficient. 

Attention  here  is  to  be  focused  primarily  on  the  afterbody,  adverse  pressure 
gradient  region.  As  illustrated  in  Figure  2,  these  two  vehicles  taken  together 
afford  a broad  range  of  adverse  pressure  gradient  over  their  quite  dissimilar 
afterbodies. 

All  boundary  layer  calculations  were  performed  with  axisymmetric, 
finite  difference  boundary  layer  computer  codes*  such  as  that  built  into  the 
Transition  Analysis  Program  System  (TAPS)  code,  which  is  now  operational  at 
Autonetics.  Detailed  stability  analyses  were  performed  by  Professor  F.  Wazzan 
of  UCLA.  These  were  compared  to  and  extended  by  in-house  computations 
utilizing  the  TAPS  code-  The  following  step-by-step  procedures  were 
(more  or  less)  employed  during  the  computational  phase  of  this  study  for 
each  vehicle. 


*We  would  like  to  acknowledge  and  express  our  appreciation  for  the  boundary 
layer  calculations  performed  by  Dr.  C.  Merkle  of  Dynamics  Technology,  Inc., 
Torrance,  California. 
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a.  Heat  Only 

A temperature  distribution  which  would  maintain  laminar  flow  over  the 
favorable  pressure  gradient,  forebody  region  was  determined  by  trial  and  error. 
Both  shape  factor  H and  detailed  stability  calculations  were  utilized  to  insure 
that  stable  laminar  flow  was  obtained. 

b.  Suction  Only 

Weak  suction  was  applied  over  unheated  forebodies  to  investigate 
the  effects  of  suction  relative  to  heating  in  a favorable  pressure  gradient. 

Shape  factor  and  stability  calculations  were  once  again  monitored. 

2.  Afterbody 

The  heating  selected  in  1 (above)  was  applied  on  the  forebody  for  all 
subsequent  afterbody  configurations  and  calculations. 

a.  Heat  Only 

Various  temperature  distributions,  up  to  a practical  limit  of 
overheat,  were  tried  in  the  afterbody  region  to  investigate  the  effects  of 
heating. 

b.  Suction  Only 

Suction  was  applied  to  the  afterbody  region  while  monitoring  the 
shape  factor  to  determine  the  minimum  levels  and  distributions  required  to 
stabilize  the  boundary  layer  to  the  90-95%  of  length  position. 
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c.  Suction  Plus  Heat 

The  best  suction  distribution  from  (b)  was  then  modified  by  removing 
circumferential  strips  of  suction  and  replacing  them  with  strips  of  constant 
temperature,  heated  surface.  This  is  shown  schematically  on  Figure  6.  Width  of 
the  heated  strips,  abruptness  of  transition  from  "suction  on"  to  "suction  off"  (and 
vice  versa)  and  the  locations  of  the  heated  strips  were  all  varied.  Initially, 
the  results  were  scrutinized  by  monitoring  the  evaluation  of  H.  Later  on, 
however,  detailed  stability  calculations  were  performed  for  the  most  interesting 
cases. 
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II.  RESULTS 

The  forebody  section  is  taken  to  be  the  forward  position  of  the  vehicle 
extending  from  the  stagnation  point  to  the  minimum  pressure  point.  Minimum 
pressure  occurs  downstream  of  maximum  diameter.  While  attention  was  focused 
on  the  afterbody  region,  calculations  were  performed  and  results  will  be 
presented  for  both  regions  of  both  vehicles.  Input  data  consists  of  geometry 
or  coordinates,  potential  flow  velocity  field,  surface  temperature  and/or 
suction  velocity  distribution.  In  raw  form,  the  results  include;  displacement 
thickness,  6*;  shape  factor,  H;  and  friction  coefficient,  Cf,  at  each  streamline 
computing  station  and  total  wake  drag  (obtained  by  a Squire-Young  integration 
through  the  downstream  boundary  layer).  Detailed  stability  characteristics  are 
presented  for  the  most  interesting  cases. 

A.  AR-3  VEHICLE 

The  hull  shape  shown  in  Figure  2a  and  designated  AR-3  was  chosen  for 
this  study  because  it  has  a large  favorable  pressure  gradient  region  which 
extends  to  approximately  72%  of  body  length.  As  such,  it  is  well  suited  to 
a laminar  flow  application.  Calculations  indicate  that  the  forebody  of 
this  hull  will  support  laminar  flow  with  a minimal  heat  input.  However, 
beyond  the  minimum  pressure  location,  the  pressure  gradient  becomes 
increasingly  adverse,  and  laminar  separation  and  subsequent  transition  to 
turbulent  flow  occur  rapidly. 


1 . Forebody  Region  (Favorable  Pressure  Gradient) 

Boundary  layer  calculations  over  the  forebody  region  of  the  AR-3  hull  have 
been  conducted  to  determine  the  minimum  temperature  distribution  which  will 
maintain  stabilized  laminar  flow.  The  process  is  basically  a trial  and  error 
approach  in  which  the  boundary  layer  shape  factor  is  monitored.  Once  a good 
profile  is  obtained,  more  detailed  and  complex  stability  calculations  are 
employed  to  verify  that  the  flow  remains  stable. 


The  final  iteration  of  temperature  distribution  for  the  AR-3  forebody  is 
presented  in  Table  1.  Also  included  in  this  table  are  geometry,  potential  flow 
and  surface  temperature  distributions  for  the  entire  body;  and  boundary  layer 
calculation  outputs  including  displacement  thickness,  6*,  shape  factor,  H,  and 
heat  transfer  coefficient,  h,  to  the  point  of  laminar  separation.  Figure  7 
shows  a plot  of  the  shape  factor,  H,  over  the  favorable  pressure  gradient  region. 
For  this  case,  labeled  "no  suction",  the  shape  factor  falls  monotonically  with 
S/L  over  the  forebody  and  then  rises  rapidly  as  the  adverse  pressure  gradient 
region  is  encountered,  S/L  > .7.  Stability  calculations  for  this  case  by 
Professor  Wazzan  of  UCLA  show  a maximum  Tollmien-Schl ichting  wave  disturbance 
amplification  of  e2,1*  in  the  favorable  pressure  gradient  region.  Thus,  the 
boundary  layer  is  expected  to  remain  laminar  and  quite  stable  to  the  point  of 
separation  at  the  onset  of  the  adverse  pressure  gradient  region. 

2 . Afterbody  Region  (Adverse  Pressure  Gradient) 


Having  defined  the  forebody  conditions  sufficient  to  maintain  laminar 
boundary  layer  flow  to  the  point  of  minimum  pressure,  subsequent  boundary  layer 
calculations  were  made  to  determine  the  extent  to  which  stabilized  laminar  flow 
could  be  carried  into  the  adverse  pressure  gradient  region.  Three  basic 
approaches  were  tried  - (a)  heat  alone,  (b)  suction  alone,  and  (c)  suction  plus 
heat.  Each  is  discussed  below. 


a.  Heat  Alone 

Extension  of  the  maximum  forebody  overheat  of  26°F  into  the  adverse 
pressure  gradient  region  results  in  the  shape  factor  history  labelled  "no  suction" 
in  Figure  7.  Greater  overheats  of  80°F,  100°F,  120°F  and  132Lf  beginning  from 
5 feet  to  13  feet  before  the  point  of  minimum  pressure  show  negligible  effect  in 
delaying  the  onset  of  laminar  separation  as  indicated  by  a zero  skin  friction. 
Figure  8 presents  a summary  of  these  cases  in  terms  of  the  skin  friction 
coefficient,  Cf,  as  a function  of  longitudinal  position.  This  plot  shows 
skin  friction  going  to  zero  and  subsequent  separation  at  approximately  the 
same  location  for  all  of  the  specified  afterbody  overheats. 
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b.  Suction  Alone 

Subsequent  boundary  layer  characteristics  were  computed  with  emphasis 
on  defining  the  level  of  afterbody  suction  which  Is  required  to  keep  the  boundary 
layer  attached  in  the  adverse  pressure-gradient  region  of  the  body.  Of  the 
several  suction  distributions  which  were  tried,  three  representative  cases  are 
given  on  Figure  9.  For  these  three  cases,  the  suction  was  commenced  at  about 
the  0.72  S/L  location  on  the  body  (S  = 39.6  feet  on  the  55-foot  body).  Beyond 
this  point,  the  suction  was  linearly  increased  until  the  .95  S/L  point,  and  then 
allowed  to  decrease  again.  "Optimization"  of  the  amount  of  suction  was  primarily 
achieved  by  rotating  the  Vw/Uo°  vs.  S/L  suction  schedule  about  the  peak  point  on 
the  curve,  although  some  checks  as  to  the  effect  of  changes  in  location  of  the 
initial  starting  point  for  the  suction  were  also  made.  With  the  "low"  suction 
schedule,  the  shape  factor  begins  to  diminish  rapidly  as  soon  as  the  suction 
section  is  encountered,  but  the  strong  unfavorable  pressure  gradient  soon 
overcomes  the  effect  of  the  suction  and  the  shape  factor  increases  very  rapidly 
until  separation  occurs.  A plot  of  the  shape  factor  versus  S/L  is  given  in 
Figure  7,  and  shows  these  effects  clearly. 

The  "intermediate"  suction  schedule  shown  on  Figure  9 kept  the 
boundary  layer  attached  slightly  longer  than  did  the  "low"  value,  but  again  as 
the  pressure  gradient  grew  more  adverse  the  boundary  layer  separated. 

The  "high"  suction  level  shown  on  Figure  9 was  sufficiently  strong  to 
keep  the  boundary  layer  attached  beyond  the  peak  in  the  suction  schedule  (S/L  = 
0.95).  The  corresponding  shape  factor  history  for  this  case  is  also  shown  on 
Figure  7.  Although  the  scale  of  the  ordinate  on  this  figure  is  highly  magnified, 
it  is  important  to  note  how  rapidly  the  integral  characteristics  of  the  boundary 
layer,  represented  by  H,  change  in  response  to  relatively  small  changes  in  the 
level  of  suction.  This  is  the  result  of  two  very  strong  competing  factors,  the 
pressure  gradient  and  the  suction.  Additional  details  of  the  behavior  of  the 
shape  factor  versus  distance  after  the  point  at  which  the  suction  is  begun  are 
given  on  Figure  10,  which  has  an  expanded  longitudinal  scale  as  compared  to 
Figure  7. 
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c.  Suction  Plus  Heat 

In  addition  to  computations  with  continuously  distributed  suction, 
some  discontinuous  suction  schedules  were  studied  to  assess  the  effectiveness 
of  surface  heating  as  a means  for  delaying  boundary  layer  separation  in 
conjunction  with  suction.  Starting  from  the  "high"  suction  schedule  of  Figure  9, 
the  initial  2.5  feet  of  suction  were  retained,  followed  by  a heated,  impermeable 
wall  (zero-suction)  region.  The  initial  2.5  feet  of  suction  was  used  to  ensure 
that  a "strong"  velocity  profile  which  was  far  from  separation  existed  at  the 
beginning  of  the  heated  non-porous  section,  and  the  length  of  the  impermeable 
section  was  varied  to  determine  how  long  the  boundary  layer  would  remain  attached 
without  suction.  Some  representative  results  for  which  the  impermeable  wall 
section  was  0.3,  0.2  and  0.1  foot  respectively  (followed  by  a return  to  the 
original  suction  schedule),  are  shown  on  Figure  11.  As  can  be  seen,  the  initial 
length  of  porous  wall  was  indeed  sufficient  to  develop  a strong  boundary  layer 
profile  which  was  far  from  separation  (H  * 1.95).  As  soon  as  the  boundary  layer 
reached  the  heated  (solid)  wall  section,  the  shape  factor  increased  dramatically 
and  separation  was  encountered  in  a distance  of  less  than  0.4  feet.  It  is 
interesting  to  note  that  after  the  suction  has  been  restored  there  is  still  a 
substantial  residual  effect  of  the  solid  wall  section  as  indicated  by  an  increase 
in  H over  that  before  the  solid  wall  section. 

\ 
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The  results  of  using  multiple  heated  sections  are  given  in  Figure  12. 

The  results  in  this  figure  are  for  an  initial  2.5  feet  of  suction  followed  by  0.1 
feet  of  heating  (solid-wall),  and  then  by  two  patterns  of  3 feet  of  suction  with 
0.2  feet  of  heating.  As  can  be  seen,  the  second  heating  pulse  showed  a peak 
shape  factor  of  nearly  2.5  while  the  third  heating  pulse  led  to  separation.  By 
comparison,  it  is  noted  that  an  initial  pulse-width  of  0.2  feet  followed  by  a 
second  pulse  of  either  0.2  or  0.1  feet  led  to  separation  at  the  second  (as 
opposed  to  the  third)  solid  wall  section. 

In  all  these  discontinuous  suction  studies,  a distance  of  0.1  foot  was 
allowed  for  the  suction  to  go  from  "on"  to  "off"  or  from  "off"  to  "on".  Thus, 
in  cases  where  the  heated  section  was  0.1  feet  long,  the  suction  was  less  than 
its  corresponding  value  for  the  continuous  suction  case  for  a distance  of  0.3  feet. 
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Finally,  calculations  for  cases  wherein  the  solid  portion  of  the 
wall  was  left  unheated  showed  almost  no  difference  from  these  heated  cases;  for 
this  body  then,  apparently  the  adverse  pressure  gradient  is  so  strong  that  it 
completely  dominates  any  favorable  effects  of  heating. 

To  enhance  the  likelihood  of  stabilized  laminar  flow,  a revised 
suction  plus  heat  distribution  incorporating  two  0.1  foot  bands  of  no  suction 
with  30°F  surface  overheat  was  employed  as  shown  in  Figure  13.  Resulting  boundary 
layer  calculations  for  this  case  are  presented  in  Figure  14  which  depicts  shape 
factor,  H,  and  skin  friction  coefficient,  Cf,  as  function  of  the  non-dimensional ized 
arc  length,  S/L.  Results  of  detailed  stability  calculations  for  this  case  by 
Professor  F.  Wazzan  of  UCLA  are  given  in  Table  2 and  show  no  disturbances 
amplifying  to  e' 5 or  greater  up  to  S/L  = .90. 

In  other  words,  in  spite  of  strong  trends  to  instability  and  transition 
(as  evinced  by  a rapidly  increasing  shape  factor)  in  regions  where  suction  was 
switched  off,  the  afterbody  flow  apparently  remained  laminar  for  the  intermittent 
suction  distribution  of  Figure  13.  The  effects  of  surface  heating  in  regions  of 
strong  adverse  gradient  (as  AR-3)  are  negligible.  The  distances  along  the  AR-3 
afterbody  over  which  suction  may  be  switched  off  and  stability  still  be 
maintained  is  dependent  only  upon  the  pressure  gradient  and  the  suction  rate 
histories.  These  distances  are  short  and  represent  extremely  small  reductions  in 
the  overall  afterbody  suction  flux. 

B.  REICHARDT  BODY 

A second  hull  with  contrasting  characteristics  to  those  of  the  AH-8  was  also 
chosen  for  this  study.  This  is  a Reichardt  body  of  5 to  1 fineness  ratio  (shown 
in  Figure  2b).  While  the  Reichardt  shape  is  more  practical  from  the  point  of  view 
of  enclosed  packaging  volume,  it  has  a pressure  gradient  distribution  which  is 
less  conducive  to  laminar  flow  stabilization  (Figure  3).  The  pressure  is  seen 
to  be  nearly  constant  over  the  middle  70  percent  of  its  length.  Moreover,  its 
pressure  gradient  becomes  adverse,  albeit  weakly,  at  approximately  35  percent 
of  length. 
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1 . Forebody  Region 

Boundary  layer  calculations  were  performed  over  the  forebody  region  of 
the  Reichardt  hull  to  determine  a temperature  distribution  which  would  maintain 
stabilized  laminar  flow  to  the  point  of  minimum  pressure.  The  result  of  several 
iterations  was  a temperature  distribution  which  rose  from  ambient  temperature 
at  the  nose  to  40°F  overheat  at  minimum  pressure  (S/L  « .39).  Details  of  the 
input  data  and  output  calculations  for  the  final  temperature  distribution  are 
given  in  Table  3.  A plot  of  shape  factor,  H,  versus  arc  length  position,  S/L, 
for  this  case,  labelled  "heated  - No  suction",  is  shown  in  Figure  15.  This 
temperature  distribution  maintains  the  boundary  layer  shape  factor  at  a nearly 
constant  value  of  2.34  over  nearly  75  percent  of  the  body,  including  nearly  65 
percent  of  the  adverse  pressure  gradient  region.  Stability  calculations  for 
this  case  show  all  disturbances  damped  to  the  minimum  pressure  location.  For 
purposes  of  comparison,  plots  of  shape  factor  for  cases  with  (a)  no  heat  or 
suction,  and  (b)  no  heat  but  small  amounts  of  suction  are  shown  in  the  same 
Figure.  In  the  no  heat  or  suction  case,  the  shape  factor  passes  through  2.4  at 
0.05  percent  of  body  length  and  remains  at  a value  of  2.6  to  70  percent  of  body 
length.  Boundary  layer  transition  is  expected  when  H grows  larger  than  about 
2.45.  The  low  level  of  suction  maintains  the  shape  factor  at  a safe,  low  value 
of  2.25  to  80  percent  of  body  length. 

2.  Afterbody  Region 

Similar  to  the  approach  outlined  above  for  the  AR-3  hull,  three  different 
systems:  (a)  heat  alone,  (b)  suction  alone,  and  (c)  suction  plus  heat  were 

employed  in  an  attempt  to  extend  stabilized  laminar  flow  from  the  minimum  pressure 
point  into  the  adverse  pressure  gradient  region.  These  are  discussed  in  turn 
below. 


a.  Heat  Alone 

Extension  of  a 40°F  overheat  from  the  minimum  pressure  location  into 
the  adverse  pressure  gradient  results  in  shape  factor,  H,  and  skin  friction 
coefficient,  Cf,  as  shown  in  Figure  16.  Beyond  an  arc  length  position  of  S/L  = 0.85 


-16- 


the  skin  friction  begins  to  drop  rapidly  and  separation  occurs  between 
S/L  = .89  and  S/L  = .92.  Boundary  layer  calculations  were  also  run  for 
different  overheats,  i.e.  0,  20,  60  and  80°F.  Figure  17  presents  the  arc  length 
variation  of  the  skin  friction  coefficient  for  all  five  cases,  and  shows  that 
differences  in  overheat  level  have  negligible  effect  as  the  pressure  gradient 
becomes  increasingly  adverse. 

b.  Suction  Alone 

In  order  to  delay  separation,  a number  of  suction  distributions  were 
applied  to  the  extreme  adverse  pressure  gradient  region  of  the  Reichardt  body. 
Boundary  layer  calculations  were  carried  out  in  an  effort  to  define  an  "optimum" 
suction  distribution  which  would  delay  separation  to  beyond  S/L  = 0.90. 

The  shape  factor  for  the  initial  region  of  the  body  was  given  in 
Figure  16.  Under  the  influence  of  the  nearly  constant  pressure  gradient,  the 
shape  factor  remains  nearly  constant  at  about  2.35  over  the  first  45  feet  of  the 
vehicle  (S/L  = 0.80).  Beyond  this  point,  the  pressure  gradient  becomes  more 
unfavorable  and  the  shape  factor  starts  to  increase.  Details  of  the  shape  factors 
on  this  aft  portion  of  the  vehicle  for  the  various  suction  distributions 
investigated  are  given  in  Figure  18.  In  the  absence  of  suction,  curve  0,  the 
shape  factor  heads  towards  separation,  but  not  so  rapidly  as  for  the  AR-3  vehicle. 
The  suction  distributions  corresponding  to  the  shape  factor  histories  shown  in 
Figure  18  are  illustrated  in  Figure  19.  When  the  maximum  suction  is  applied,  the 
shape  factor  drops  rapidly  (curve  1).  Note  that  the  peak  in  this  suction  curve  is 
the  same  as  the  one  used  for  the  AR-3  vehicle,  but  that  the  axial  extent  of  the 
suction  is  considerably  shorter.  Successive  reductions  in  the  suction 
distribution  are  given  by  Schedules  2,  3,  and  4 on  Figure  19,  with  the 
corresponding  changes  in  -the  shape  factor  distribution  given  on  Figure  18.  Each 
of  these  suction  distributions  is  sufficient  to  ensure  that  the  boundary  layer 
remains  attached,  and  the  latter.  Schedule  4,  has  been  chosen  as  the  "optimum" 
distribution  since  further  reductions  in  suction  lead  to  early  separation. 
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c.  Suction  Plus  Heat 


As  in  the  AR-3  calculations,  this  "optimum"  continuous  suction 
distribution  was  them  made  discontinuous  by  interspersing  impermeable,  heated 
surface  sections  to  determine  the  feasibility  of  using  surface  heating  to  keep 
the  boundary  layer  attached.  Because  the  adverse  pressure  gradient  on  this 
Reichardt  body  is  not  as  steep  as  that  on  the  AR-3  vehicle,  the  lengths  of  the 
impermeable  sections  could  be  made  considerably  longer.  The  shape  factor 
distributions  for  the  discontinuous  suction  distributions  are  also  given  in 
Figure  18.  Curve  number  5 corresponds  to  an  initial  run  of  2 feet  of  porous 
surface  followed  by  an  impermeable  heated  surface  over  the  rest  of  the  vehicle. 
(The  heating  used  for  the  Reichardt  body  was  AT  = 40°F).  Following  the  two  feet 
of  suction,  the  heated  surface  keeps  the  boundary  layer  attached  for  about  3 feet 
before  separation  occurs.  In  curve  number  6 (on  Figure  19),  suction  was  applied 
in  two  two-foot  intervals  separated  by  a 0.4  foot  heated  section  with  the  entire 
rear  surface  of  the  vehicle  again  being  heated.  Here,  separation  occurred 
rapidly  once  the  final  heated  surface  was  reached. 

Finally,  the  results  for  a discontinuous  suction  distribution  with 
two  0.4  foot  intervals  of  heating  is  given  on  Figure  20  (labeled  "hot").  The 
suction  distribution  for  this  case  is  also  given  on  Figure  19.  This  last  case 
follows  the  suction  levels  of  Schedule  4 with  the  inclusion  of  the  two  "suction 
off"  bands  of  Schedule  6.  In  this  case,  the  shape  factor  increased  slightly  over 
the  heated  surfaces,  but  the  boundary  layer  remained  attached.  The  curve  labeled 
"cold"  on  Figure  20  corresponds  to  tha  same  suction  distribution  but  with  the 
heat  turned  off  in  the  solid  wall  regions.  For  comparison,  the  continuous 
suction  distribution  is  also  given  (lowest  curve).  This  comparison  shows  that 
the  heating  does  provide  a measure  of  protection  against  separation  on  the 
impermeable  portions  of  the  surface  for  these  more  moderate  adverse  pressure 
gradients  (compared  to  the  AR-3),  but  that  after  a distance  the  history  effects 
die  out. 
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The  suction  plus  heat  scheme  which  gave  the  shape  factor  history  of 
Figure  20  was  analyzed  for  its  corresponding  boundary  layer  stability 
characteristics.  The  shape  factor  and  friction  coefficient  for  this  case, 
designated  as  1-W,  have  been  replotted  in  Figure  21.  The  calculated  stability 
growth  factors  are  presented  in  Figure  22.  This  is  a plot  of  the  natural 
logarithm  of  the  predicted  disturbance  amplitude  ratios  for  each  critical  frequency 
invesitgated  as  a function  of  arc  length.  The  term  growth  factor  is  defined  as. 


and  the  frequency  has  been  non-dimensional i zed  in  the  usual  way;  i.e. 

* wv 
w*  = — 

U2 

where  w is  the  dimensional  Tollmien-Schl ichting  wave  frequency,  v is  the  kinematic 
viscosity  and  U is  the  vehicle  speed. 

The  growth  factor  is  seen  in  Figure  22  to  increase  with  arc  length 
to  a maximum  for  each  of  the  critical  frequencies.  Progressing  further  along  the 
body  beyond  the  position  of  maximum  growth  factor,  the  waves  are  expected  to  damp 
and  the  growth  factor  to  decrease.  This  aspect  is  not  investigated  by  the  computer 
boundary  layer  stability  code,  nor  is  it  shown  on  the  figure.  These  results 
indicate  a growth  factor  of  nearly  8 (for  the  non-dimensional  frequency  .125  x 10* s) 
which  is  considered  large  enough  to  indicate  a strong  possibility  of  transition. 

A series  of  modifications  to  this  suction  plus  heat  distribution  were 
then  made  in  an  attempt  to  decrease  the  maximum  growth  factor  and  the  probability 
of  transition.  A description  of  these  modifications  is  presented  in  Table  4 and 
will  be  described  here  along  with  their  rationale. 


i 
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Case  2-W  represents  a modification  of  Case  1-W  in  which  heat  was 
added  to  the  porous,  suction  walls  of  Case  1-W.  Specifically,  a 40°F  wall 
overheat  was  applied  in  the  regions  of  suction.  Shape  factor  and  skin  friction 
coefficient  histories  are  shown  in  Figure  23;  corresponding  growth  factor 
histories  are  presented  in  Figure  24.  Comparison  to  Case  1-W,  Figures  21  and  22, 
reveals  only  slight-to-negligible  differences  between  the  two  cases  due  to  the 
addition  of  heat. 

The  most  amplified  disturbance  frequencies  for  both  Case  1-W  and 
Case  2-W  are  seen  to  grow  to  unacceptibly  large,  unconservative  levels  approaching 
the  "e9"  transition  prediction  level.  In  an  effort  to  reduce  these  levels  and 
improve  the  stability  of  the  attendant  boundary  layer,  a low  level  of  suction 
was  applied  upstream  of  that  for  Cases  1-W  and  2-W.  The  results  for  this 
modification,  Case  3-W,  are  shown  in  Figures  25  and  26.  For  this  configuration, 
disturbances  at  dimensional  frequencies  in  the  vicinity  of  0.1-0.15  x 10' 5 are 
seen  to  grow  well  above  e9  near  the  arc  length  location  at  which  suction  commences. 

This  is  a rather  curious  result  in  that  suction  which  is  generally  thought  to  be 
stabilizing  is  seen  here  to  be  destabilizing.  The  explanation  for  this  observation 
lies  in  the  often  overlooked  fact  that  suction,  which  both  thins  and  distorts  the 
boundary  layer  profile,  may  increase  a boundary  layer  sensitivity  to  a disturbance 
wave.  The  characteristics  of  the  most  amplified  disturbances  are  established  by  the 
thickness  and  shape  of  the  boundary  layer  as  well  as  by  the  freestream  velocity. 

Suction,  therefore,  has  the  power  to  create  a more  unstable  environment  for 
disturbances  propagating  from  some  upstream  location  by  essentially  "tuning"  the 
boundary  layer  to  waves,  which,  without  suction,  would  soon  decay.  This 
phenomenon  is  displayed  in  Figure  26  by  the  growth  factor  history  for  the 
disturbance  wave  w*  = 0.15  x 10~6  which  begins  to  flatten  out  at  S/L  = .68  and 
then  amplifies  very  rapidly  near  S/L  = .70  where  suction  is  turned  on. 
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For  Case  4-W,  the  upstream  suction  levels  were  increased  over  those 
of  Case  3-W.  Figures  27  and  28  show  lower  shape  factor,  higher  skin  friction 
coefficient  and  reduced  growth  factors  compared  to  Cases  1-W,  2-W  and  3-W. 

Case  4-W  shows  particularly  good  stability  characteristics  with  a maximum 
disturbance  growth  factor  of  less  than  4.  All  disturbances  are  rapidly  and 
effectively  damped  (not  well  displayed  in  the  figure)  in  the  presence  of  this 
suction  distribution. 

For  the  fifth  and  final  iteration,  the  upstream  suction  levels  were 
reduced  to  values  intermediate  to  those  of  Cases  3-W  and  4-W  (reference  Table  4). 
The  results,  displayed  in  Figures  29  and  30,  are  extremely  interesting, 
particularly  in  regard  to  the  previous  comments  on  boundary  layer  "tuning"  by 
suction.  Figure  30,  which  presents  the  growth  factor  histories  for  the  most 
amplified  disturbance  frequencies,  gives  evidence  of  suction's  ability  both  to 
tune  and  to  detune  a boundary  layer  to  a given  disturbance  wave.  The  disturbance 
w*  = 0.17  x 10" 5,  for  example,  is  seen  to  be  apparently  approaching  its  maximum 
growth  at  S/L  = .68  only  to  experience  a renewed,  very  rapid  growth  at  the 
commencement  of  suction  near  S/L  = .7.  For  S/L  £ .7,  however,  the  same  wave 
suddenly  begins  to  decay  very  rapidly  due  to  detuning  of  the  boundary  layer 
by  the  specified  suction  distributions. 

These  results  are  re-examined  in  the  light  provided  by  existing 
analytical  studies  and  discussed  at  greater  depth  in  Section  IV. 


III.  HULL  DRAG  ANALYSIS 
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To  provide  a measure  of  the  effectiveness  of  the  different  laminar  flow 
stabilization  schemes  investigated,  calculations  of  total  hull  drag  were  carried 
out.  That  drag  should  be  the  basis  for  this  comparison  instead  of,  say,  cost 
or  overall  system  reliability  or  simplicity  is  a moot  point  at  the  present  time. 
Ultimately  all  relevant  factors  would  have  to  be  taken  into  account.  However, 
estimates  of  drag  are  relatively  simple  to  perform  and  do  provide  information 
on  vehicle  power  requirements,  which  is  fundamental.  The  results  of  these 
calculations,  which  are  summarized  below,  enable  a first  order  comparison  not 
only  among  the  different  stabilization  schemes,  but  also  between  the  two  hull 
shapes. 

Several  simplifying  assumptions  have  been  made.  In  particular,  the  vehicles 
are  assumed  to  be  without  appendages  such  as  fins  or  control  surfaces  so  that 
only  hull  drag  calculations  need  be  performed.  The  total  hull  drag  is  assumed 
to  be  comprised  of  (a)  wake  drag,  (b)  drag  equivalent  of  suction  power,  and 
(c)  drag  equivalent  of  heating  power.  Each  drag  component  is  expressed  as  a 
drag  coefficient  in  the  conventional  manner,  defined  below.  The  total  hull  drag 
coefficient,  is  expressed  as. 


where  Cpw  is  the  wake  drag  coefficient,  C[js  is  the  equivalent  suction  drag 
coefficient,  and  Cph  is  the  equivalent  heating  drag  coefficient. 


The  wake  drag  contribution  was  calculated  from  the  momentum  loss  integral 
at  the  final  hull  station,  as  given  by  the  boundary  layer  computer  code,  according 
to  the  Squire- Young  formula 


Here  U is  the  local  velocity,  U«>,  is  the  freestream  value,  A is  the  reference 
area  - generally  taken  to  be  the  wetted  hull  area  - and  S0  is  the  toroidal  area 
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bounded  by  the  hull  radius  and  the  hull  radius  plus  the  boundary  layer 
momentum  thickness.  A correction  was  applied  for  the  non-ambient  static 
pressure  acting  at  the  final  hull  station. 


The  equivalent  suction  drag  coefficient,  C[)s,  was  determined  from  the  power 
requirements  of  the  suction  pump.  Assuming  the  case  of  a single  stage  suction 
pump,  the  pump  power,  Ps,  may  be  expressed  as. 


P$  = Q • Ap 

where  Q is  the  total  volume  flow  rate  of  sucked  fluid  and  Ap  represents  the 

total  pressure  rise  through  the  suction  system.  The  suction  flow  rate,  Q,  is 

obtained  by  integrating  the  suction  velocity  over  the  suction  area.  From  the 

single  stage  pump  assumption,  the  pump  must  be  sized  to  provide  pressure  rise 

from  the  lowest  local  static  pressure  value  (below  ambient)  on  the  -jction 

surface  to  the  freestream  dynamic  pressure,  qa>,  (above  ambient)  while  additionally 

overcoming  pressure  drops  through  the  suction  skin.  Experience  has  shown  that 

for  operational  systems,  skin  pressure  drops  are  generally  on  the  order  of  3%  of 
1 fi 

dynamic  pressure.  Ducting  is  assumed  large  enough  to  produce  negligible  pressure 
losses  and  pump  efficiency  (electrical  to  mechanical  conversion)  is  assumed  ideal. 
The  suction  pump  power,  Ps,  is  converted  to  a non-dimensional  equivalent  drag 
coefficient  by  dividing  by  the  product  of  freestream  velocity,  dynamic  pressure 
and  wetted  hull  area;  as, 


P 

s 


Uco  (H  p[)i) 


A 


The  equivalent  heating  drag  coefficient,  was  determined  in  a similar 
manner.  Waste  heat  from  the  vehicle  propulsion  system  was  assumed  to  be  the 
source  of  heating.  Therefore,  the  only  penalty  charged  against  equivalent  drag 
due  to  heating  results  from  the  transport  and  distribution  system.  A configuration 
suitable  for  transporting  and  distributing  heat  to  the  vehicle  surface  was 
proposed  in  Reference  15.  It  consists  of  metallic  tubes  (copper  or  nickel)  which 
are  embedded  in  a fiberglass  vehicle  inner  hull  and  electroformed  onto  an 
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outer  metallic  skin  and  through  which  heating  fluid  is  pumped.  The  pump 
power,  Ph,  required  to  distribute  waste  heat  via  this  secondary  fluid  - heat 
exchanger  scheme  may  be  estimated  as  the  product  of  secondary  fluid  volume 
flow  rate  and  a fixed  pressure  loss. 


This  concept  of  surface  heating  has  been  investigated  in  detail  for  a 
vehicle  shape  similar  to  the  AR-3  at  Autonetics.^  The  vehicle  surface 
temperature  distributions  and  total  heat  transfer  requirements  for  the 
contour  studied  are  similar  to  those  of  both  the  AR-3  and  the  (5:1) 
Reichardt  vehicles  under  study  here.  Estimates  of  pump  power  by  means  of 
required  flow  quantity  and  ducting  pressure  losses  may,  therefore,  be 
obtained  from  the  results  of  the  study  of  Reference  17.  Heating  fluid  flow 
rates  on  the  order  of  180  Lbm/sec  (or  2.8  ft3/sec)  and  ducting  pressure 
losses  on  the  order  of  5 psi  (or  720  psf)  have  been  shown  to  be  adequate. 
Pressure  losses  were  subsequently  increased  by  50  percent  to  account  for 
manifolds  and  connections.  The  flow  rate,  pressure  loss  product  was 
non-dimensionalized  by  the  freestream  velocity  - dynamic  pressure  - wetted 
hull  area  product  to  generate  the  equivalent  heating  drag  coefficient;  as, 

C . li 

°h  {«,  p u')  A 
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This  coefficient  was  divided  into  two  components  relating  to  the  forebody  and 
afterbody  heating  systems.  The  afterbody  contributions  will  be  seen,  however, 
to  be  negligible  in  comparison  to  all  others. 


Following  these  procedures,  total  drag  estimates  were  made  for  the  six  most 
interesting  cases  where  heat-plus-suction  was  applied  over  the  vehicle  afterbody. 
The  results  are  summarized  in  the  following  table. 
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CONFIGURATION 

CD 

w 

C°s 

C°h 

(FOREBODY) 

c“h 

(AFTERBODY) 

AR-3 

2.7  x 10"5 

9.0  x 10"5 

.35  x 10"5 

1.1  x 10'10 

12.1  x 10‘5 

Reichardt  (1) 

6.3  x 10'5 

1.5  x 10"5 

.42  x 10'5 

1.3  x 10"9 

8.2  x 10‘5 

Reichardt  (2) 

6.3  x 10~5 

1.5  x 10'5 

.42  x 10"5 

1.3  x 10"9 

8.2  x 10"5 

Reichardt  (3) 

5.5  x 10'5 

2.4  x 10'5 

.28  x 10"5 

1.3  x 10"9 

8.2  x 10~5 

Reichardt  (4) 

4.2  x IQ'5 

4.2  x 10'5 

.28  x 10'5 

1.3  x 10"9 

8.7  x 10"5 

Reichardt  (5) 

5.2  x 10-5 

2.8  x 10'5 

.28  x 10’5 

1.3  x ID'9 

8.2  x 10'5 

Under  the  assumption  that  there  is  no  drag  penalty  due  to  heat  generation, 
the  total  equivalent  heating  drag  coefficient  is  apparently  a very  small  fraction 
of  the  total  hull  drag  coefficient.  Moreover,  the  afterbody  contribution  to  the 
total  equivalent  heating  drag  coefficient  is  at  least  five  orders  of  magnitude 
smaller  than  the  other  contributions  and  is  therefore  negligible.  A further 
observation  is  that  for  the  Reichardt  body,  the  variations  in  afterbody  suction 
distribution  which  were  investigated  produced  only  a 5 percent  change  in  the  total 
hull  drag  coefficient.  This  is  due  to  the  following  phenomenon:  while  larger 

suction  coverage  increases  the  equivalent  suction  drag  coefficient,  the  wake  drag 
coefficient  is  simultaneously  decreased  as  the  boundary  layer  is  sucked  off  to 
diminish  the  net  effect.  Previous  experience  with  suction  systems  for  boundary 
layer  control  has  demonstrated  this  phenomenon:  for  low  suction,  the  latter 

effect  dominates  the  former  and  Cpt  decreases  as  suction  rate  increases.  For 
larger  suction,  the  reverse  is  true;  that  is,  Cpt  increases  as  suction  rate 
increases.  Further  discussion  on  this  interesting  point  is  contained  in  the 
following  section. 
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IV.  DISCUSSIONS  AND  CONCLUSIONS 


The  effects  of  heating  and  suction  may  be  conveniently  studied  on  the 
Re^*  - H plane,  which  gives  a useful  but  approximate  picture  of  the  stability 
characteristics  of  the  boundary  layer.  Results  for  the  afterbody  of  the  AR-3 
body  (Reference  Figure  14)  are  plotted  on  Figure  31  along  with  the  critical 
Reynolds  number  curve.  All  flow  states  to  the  left  of  the  critical  curve  admit 
no  growing  disturbances;  flow  states  to  the  right  are  unstable.  It  is  expected 
that  a flow  which  suffers  transition  by  the  mechanism  of  linear 
Tol lmien-Schl ichting  waves  growing  ever  larger  will  have  a trajectory  on  the 
Re^*  - H plane  which  extends  to  the  right  of  the  critical  curve.  The  trajectory 
for  the  AR-3  body  begins  at  the  minimum  pressure  point  and  proceeds  as  indicated 
by  the  arrows  through  the  suction  (to  lower  H)  and  heating  (to  higher  H)  regions. 

As  suction  is  applied,  the  boundary  layer  moves  away  from  the  critical  curve, 
but  returns  rapidly  toward  it  when  the  first  heating  pulse  (no  suction)  is 
applied.  After  the  suction  has  once  more  made  the  boundary  layer  strongly  stable, 
the  second  heating  pulse  forces  the  boundary  layer  into  the  unstable  range. 

For  the  continuous  suction  case,  the  boundary  layer  remains  well  below  the 
critical  Re^*  curve,  indicating  that  the  suction  distribution  selected  is 
sufficient  to  prevent  transition  (and  separation). 

The  Reichardt  body.  Case  1-W,  fares  considerably  worse  according  to  this  Re{*  - 
H plane  representation  (see  Figure  32).  Shortly  downstream  of  minimum  pressure, 
the  boundary  layer  becomes  unstable  and  moves  continuously  further  from  the  critical 
line.  The  applied  suction  draws  it  back  towards  the  stable  region,  and,  even  in 
this  discontinuous  case  eventually  stabilizes  the  boundary  layer.  In  order  to  keep 
this  body  laminar,  indications  are  that  the  suction  would  have  to  be  applied  much 
earlier  in  the  constant  shape  factor  region  near  the  front  of  the  body.  This  was 
demonstrated  in  Section  II  by  test  Cases  3-W,  4-W  and  5-W.  Whether  or  not  a 
reasonable  amount  of  suction  near  the  front  of  the  body  would  keep  the  boundary 
layer  stable  over  the  entire  (almost)  constant-pressure  region  is  not  clear  from 
this  presentation.  What  is  quite  clear  are  the  trends  of  a boundary  layer 
trajectory  on  this  plane  under  the  action  of  either  suction  or  heating.  Suction 
is  apparently  capable  of  shifting  the  trajectory  to  much  lower  values  of  H and 
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increased  boundary  layer  stability.  Practical  amounts  of  heating,  on  the 
other  hand,  for  the  levels  of  pressure  gradient  encountered  with  the  two  bodies 
under  study  here  are  much  less  effective  in  this  regard.  These  calculations 
and  observations  suggest  that  the  relative  effectiveness  of  the  two  mechanisms 
can  be  monitored  by  trajectory  changes  in  the  Re^*  - H plane. 

Further  insight  is  provided  by  a detailed  investigation  of  the  boundary  layer 
velocity  profile  and  its  second  derivative.  These,  in  fact,  represent  the 
boundary  layer  as  input  into  the  linear  stability  analysis.  Of  particular 
interest  is  the  evolution  of  the  velocity  and  second  derivative  profiles  in  a 
region  where  suction  is  being  turned  off  and  heating  is  being  turned  on.  Four 
stations  through  such  a region  have  been  selected  on  the  afterbody  of  the  AR-3 
body,  located  by  number  on  Figure  12.  Velocity  profiles  for  Station  1,  which 
is  in  a suction-only  region,  and  for  Station  3,  which  is  in  a heat-only  region 
are  displayed  in  Figure  33.  The  former  is  seen  to  be  fuller  and  thinner  than  the 
latter,  which  reflects  both  lower  H and  increased  stability  for  the  suction  over 
the  heated  boundary  layer.  Thinner,  more  energetic  boundary  layers  are  at  once 
more  stable  and  better  able  to  negotiate  adverse  pressure  gradients. 

The  second  derivative  profiles  for  the  four  stations  are  plotted  on  Figure  34. 
As  suction  is  switched  off  and  heating  is  switched  on,  the  second  derivative 
profile  is  seen  to  suffer  a zero  crossing  just  before  Station  3.  Rayleigh's 
Theory  and  experience  have  demonstrated  that  a zero  crossing  in  the  second 
derivative  profile  (which  corresponds  to  an  inflection  point  in  the  velocity 
profile)  signals  the  onset  of  instability  within  the  flow.  The  boundary  layer, 
therefore,  goes  unstable  in  passing  from  the  region  of  suction  to  that  of  heating 
in  this  adverse  pressure  gradient.  The  width  of  the  heated  or  no  suction  region 
must  be  less  than  the  distance  for  instability  to  develop  into  transition. 

That  this  is  so  has  been  demonstrated  by  detailed  stability  calculations  for  the 
same  case,  which  indicate  that  transition  does  not  occur. 
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Another  presentation  which  adds  insight  into  the  mechanisms  of  suction 
and  heating  and  their  effectiveness  in  an  adverse  pressure  gradient  was 
considered  by  Aroestyand  Berger  (Reference  10).  In  an  approximate  analysis  of 
(for  Falkner-Skan  water  flows)  suction  and  heating  in  the  control  of  laminar 
separation,  they  were  able  to  develop  relationships  which  indicate  the  level 
of  heating  and/or  suction  required  to  delay  separation  in  an  adverse  pressure 
gradient.  Their  results  (their  Figures  3 and  4)  are  included  here  for  reference 
in  Figure  35.  These  plots  represent  the  minimum  amount  of  heating  or  suction 
necessary  to  maintain  an  attached  laminar  boundary  layer  in  an  adverse  pressure 
gradient.  Points  to  the  right  (below  the  curves)  represent  separated  and  to 
the  left  (above)  represent  unseparated  flows.  At  a length  Reynolds  number  of 
250  x 106,  a 40°F  surface  overheat  has  about  the  same  effectiveness  in  delaying 
laminar  separation  as  a level  of  suction  vw/Uoo  = 0.006  x 1 0" 3 ! Thus,  surface 
heating  is  predicted  to  be  of  limited  utility  and  of  greatly  reduced  effectiveness 
in  delaying  laminar  separation,  as  compared  to  suction.  The  same  general 
conclusions  can  be  drawn  concerning  heating  and  suction  and  transition. 

The  suction  parameter  vs  6 trajectory  for  the  afterbodies  of  the  AR-3 
(reference  Figure  13)  and  the  Reichardt,  Case  1-W  (reference  Figure  9)  have  been 
plotted  in  Figure  36.  Also  plotted  is  the  separation-limit  curve  of  Figure  35 
(from  Reference  11).  Ignoring  the  downward  spikes  which  correspond  to  the  heated, 
no  suction  regions,  these  plots  illustrate  that  the  suction  levels  which  delay 
transition  for  both  the  AR-3  and  Reichardt  bodies  are  consistent  with  those 
predicted  for  delay  of  separation.  It  is  tempting  to  investigate  the  existence 
of  a transition-limit  curve  in  this  plane  which  would  indicate  the  minimum  levels 
of  suction  which  would  delay  transition  for  a given  adverse  pressure  gradient. 

Such  an  investigation  is,  unfortunately,  beyond  the  scope  of  the  present  study. 

A similar  thought  led  to  the  construction  of  the  transition  Reynolds  number  curve 
in  the  Re^*  vs  H plane  (reference  Section  II  and  Figure  5). 
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Falkner-Skan  B and  skin  friction  coefficient  are  plotted  as  a function  of 
non-dimensional  arc  length  S/L  for  the  AR-3  and  Relchardt  (Case  1-W)  bodies 
on  Figures  37  and  38,  respectively.  The  downward  spikes  in  Cf  denote  the  heated, 
no  suction  regions  of  the  afterbody.  The  results  of  Figure  35  (from  Reference  11) 
indicate  that  a wall  overheat  of  90°F  would  be  required  to  keep  a laminar 
boundary  layer  attached  for  a 8 of  -0.3.  Figure  36  shows  that  the  B's  in  the 
regions  of  the  heated,  no  suction  afterbody  strips  on  the  AR-3  body  are  on  the 
order  of  -2.0  to  -5.0.  Such  a flow  is  clearly  on  its  way  to  separation  with 
negligible  help  from  the  applied  overheat  of  30°F.  The  situation  is  somewhat 
improved  for  the  Reichardt  body.  Figure  38,  where  8 in  the  heated  strips  is  seen  to 
be  on  the  order  of  -0.2  to  -1.0.  There  is,  therefore,  some  chance  of  a little 
\ help  from  surface  heating,  for  sufficiently  large  AT,  in  the  initial  region  of  the 

Reichardt  vehicle  afterbody.  These  statements  which  are  predicated  on  results 
for  laminar  separation  rather  than  for  transition  are  entirely  consistent  with  the 

previous  discussions  and  conclusions  based  on  study  of  the  boundary  layer 

K 

trajectory  in  the  Re^*  - H plane. 

In  summary,  the  use  of  surface  heating  in  the  afterbody  region  of  the  AR-3 
or  Reichardt  (5:1  fineness  ratio)  vehicles  is  clearly  from  very  minimal  to 
absolutely  zero  utility  for  transition  suppression.  There  is  some  hope  for 
heating  in  adverse  pressure  gradients  small  compared  to  those  typical  of  vehicle 
afterbody  regions.  Boundary  layer  calculations  and  stability  analyses  performed 
in  the  present  study  indicate  the  effectiveness  of  suction  for  transition  (and 
separation)  suppression  in  an  adverse  pressure  field.  These  results  are  in 
agreement  with  the  limited  (approximate)  existing  theoretical  analyses  and  with 
, intuition.  Moreover,  the  suction  field  need  not  be  continously  distributed 

in  order  to  effectively  delay  transition.  Interruptions  in  the  suction 
distribution  which  are  short  compared  to  the  distances  necessary  for.,  small 
disturbances  to  grow  sufficiently  larger  to  cause  transition  are  permissible.  The 
width  of  such  interruptions  as  well  as  the  distribution  of  suction  intensity 
before  and  after  an  interruption  must  be  carefully  analyzed. 
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The  Re^*  - H plane  representation  provides  a useful,  convenient  and 
inexpensive  but  approximate  means  of  assessing  the  effectiveness  of  an 
intermittent  suction  field.  The  ultimate  assessment,  however,  lies  in  detailed 
analysis  of  the  boundary  layer  stability  characteristics.  This  is  especially 
so  because  of  the  ability  of  suction  to  "tune"  and  "detune"  a boundary  layer  to 
specific  disturbance  frequencies  by  modifying  the  boundary  layer  thickness  and 
velocity  profile.  While  suction  is  apparently  capable  of  stabilizing  any  laminai 
boundary  layer  flow,  the  required  suction  flow  rates  and/or  flow  distributions 
may  become  impractical. 

As  the  suction  flow  rates  increase,  so  does  the  cost  in  terms  of  equivalent 
i drag  or  energy  consumption  to  provide  the  suction  flow.  This  point  is  illustrat< 

in  Figure  39,  taken  from  Reference  18  for  an  8:1  Reichardt  body  in  air,  which 
shows  wake  drag  coefficient,  equivalent  suction  drag  coefficient  and  overall  drai 
coefficient*  as  a function  of  suction  flow  rate.  The  equivalent  suction  drag  i< 
seen  to  over  take  the  gains  realized  by  means  of  suction  in  wake  drag;  the  overa’ 
drag  coefficient  which  decreases  for  lower  levels  of  suction,  increases  for  high< 
levels.  This  point  is  also  illustrated  by  the  calculations  and  results  of 
Section  III. 

I 

Comparing  the  tabulated  results  of  wake  drag,  equivalent  suction  drag  and 

i 

overall  drag  coefficients  for  the  two  vehicles,  the  overall  drag  for  the  AR-3  is 
seen  to  be  larger  than  that  for  the  Reichardt  body.  The  wake  drag  for  the  AR-3 
’ is  considerably  smaller  than  that  for  the  Reichardt  body.  The  difference  is 

clearly  due  to  the  larger  equivalent  suction  drag  for  the  AR-3.  in  other  words, 
the  higher  levels  of  suction  required  to  stabilize  the  AR-3  afterbody  as  compare! 
to  the  Reichardt  afterbody  correspond  to  a higher  drag,  inspite  of  the  correspom 
reduction  in  wake  drag.  These  factors  need  to  be  taken  into  account  by  engineer* 
and  designers  seeking  to  find  practical  applications  for  laminar  flow  vehicles. 


*These  coefficients  are  defined  in  Section  III. 

■ ' 
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The  objective  of  this  report  was  to  present  the  results  of  a study  of 
the  stabilizing  effects  of  surface  heating  and  suction  in  the  afterbody  region 
of  two  specific  vehicles.  The  utility  of  suction  in  this  regard  was  clearly 
demonstrated.  Before  further  progress  can  be  made  in  this  area,  some  little 
understood  but  imoortant  aspects  of  the  problem,  which  emerged  during  the 
present  study,  need  to  be  investigated.  These  include; 

1.  A study  of  the  relationship  between  level  of  suction,  the  Falker-Skan  6 
parameter  and  transition  leading  to  a "transition  curve"  similar  to  the 
"separation-limit  curve"  on  Figure  35. 

2.  A detailed  study  of  the  mechanism  whereby  a boundary  layer  under  the 
action  of  suction  may  be  "tuned"  or  "detuned"  to  disturbance  waves. 

and  finally, 

3.  A thorough,  parametric  study  of  the  optimization  of  suction,  for  delay 
both  of  transition  and  of  separation,  on  laminar  vehicle  afterbodies. 

Both  continuous  and  discontinuous  suction  distributions  need  to  be 
investigated  to  insure  the  most  stable  laminar  flow.  Further 
optimization  in  terms  of  total  drag  (as  described  herein)  would  lead 
to  fully  hydrodynamically  optimized  vehicle  designs. 

Experience  has  shown  that  these  will  best  be  approached  by  careful 
theoretical  analysis  guided  and  supported  by  detailed  laboratory 
experiments.  • 
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Figure  1.  The  Effects  of  Geometry,  Heat,  and  Heat  Plus 
Suction  on  Vehicle  Performance 
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Critical  and  Transitional  Reynolds  Numbers 
(Based  on  "e9")  vs.  Shape  Factor 
(from  Reference  11) 
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Schematic  of  Afterbody  Sucti 
Plus  Heating 


Figure  11.  Discontinuous  Suction  - The  Effects  of  Varying  Width 

Heating  (AR-3) 


Discontinuous  and  Continuous  Suction 
(Reference  Page  27  and  Figuri 
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Figure  17.  Skin  Friction  Coefficient  vs.  S/L  for  the  Afterbody  Region  of  the 

Reichardt  Body  for  Several  Overheats 


Figure  18.  Shape  Factor  vs.  S/L^Reichardt  Body  - Various  Suction  Distributions 


Temperature  Distributions  Investigated  on  the  Reichardt  Body 


Effects  of  Surface  Heating  on  Shape  Factor  in  Presence  of  Unfavorabl 

(Reichardt  Body) 


21.  Shape  Factor  and  Friction  Coefficient  vs.  S/L  (Reichardt  Body,  Case  1-W) 


for  Reichardt  Case  1-W 


Reichardt,  Case  2-W 
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Factors  for  Reichardt  Body,  Case  5-W 
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Figure  31.  Shape  Factor  - Res*  History  Discontinuous  Suction  Plus  Heating 
, (AR-3  Body) 
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Figure  39.  Drag  Coefficients  vs.  Suction  Quantity 
(8:1)  Reichardt  Body  in  Air 
ReL  = 18.2  x 106 
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Table  2.  Tollmien-Schl ichting  Wave  Growth  Factors  (&n  a/a0)  AR-3  Body 
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Table  3.  Reichardt  Contour  (5:1) 
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Table  4.  Suction  Plus  Heat  Cases  Investigated  for  the  Reichardt  Body 
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where  p is  density,  v is  suction  velocity  and  is  freestream  velocity. 
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